On the origin of reentrance in 2D Josephson Junction Arrays 
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A comparative study of the magnetic properties of 
shunted and unshunted two-dimensional Josephson junc- 
tion arrays (2D-JJA) is presented. Using a single- 
plaquette approximation of the 2D-JJA model, we were 
able to successfully fit all our experimental data (for the 
temperature, AC and DC field dependencies of suscepti- 
bility) and demonstrate that the dynamic reentrance of 
AC susceptibility is directly linked to the value of the 
Stewart-McCumber parameter (3c- Based on extensive 
numerical simulations, a phase diagram (3c — (3l is plot- 
ted which demarcates the border between the reentrant 
and non-reentrant behavior in the arrays. 



I. INTRODUCTION 

According to the current paradigm, paramagnetic 
Meissner effect (PME) [1-6] can be related to the pres- 
ence of 7r-junctions [7], either resulting from the pres- 
ence of magnetic impurities in the junction [8,9] or from 
unconventional pairing symmetry [10]. Other possible 
explanations of this phenomenon are based on flux trap- 
ping [11] and flux compression effects [12] including also 
an important role of the surface of the sample [3]. Be- 
sides, in the experiments with unshunted 2D-JJA, we 
have previously reported [13] that PME manifests itself 
through a dynamic reentrance (DR) of the AC magnetic 
susceptibility as a function of temperature. These results 
have been further corroborated by Nielsen et al. [14] and 
De Leo et al. [15] who argued that PME can be simply 
related to magnetic screening in multiply connected su- 
perconductors. So, the main question is: which parame- 
ters are directly responsible for the presence (or absence) 
of DR in artificially prepared arrays? 

Previously (also within the single plaquette approxi- 
mation), Barbara et al. [13] have briefly discussed the 
effects of varying (3l on the observed dynamic reentrance 
with the main emphasis on the behavior of 2D- JJA sam- 
ples with high (and fixed) values of (3c- However, to our 
knowledge, up to date no systematic study (either ex- 
perimental or theoretical) has been done on how the (3c 
value itself affects the reentrance behavior. In the present 



work, by a comparative study of the magnetic proper- 
ties of shunted and unshunted 2D- JJA, we propose an 
answer to this open question. Namely, by using experi- 
mental and theoretical results, we will demonstrate that 
only arrays with sufficiently large value of the Stewart- 
McCumber parameter [3c will exhibit the dynamic reen- 
trance behavior (and hence PME). 



II. EXPERIMENTAL RESULTS 

To measure the complex AC susceptibility in our ar- 
rays we used a high-sensitive home-made susceptometer 
based on the so-called screening method in the reflection 
configuration [16-18]. The experimental system was cal- 
ibrated by using a high-quality niobium thin film. Pre- 
viously [18], we have shown that the calibrated output 
complex voltage of the experimental setup corresponds 
to the complex AC susceptibility. 

To experimentally investigate the origin of the reen- 
trance, we have measured x'(T) for three sets of shunted 
and unshunted samples obtained from different makers 
(Westinghouse and Hypress) under the same conditions 
of the amplitude of the excitation field h ac (ImOe < 
h ac < lOOe), external magnetic field Hd c (0 < Hdc < 
500Oe) parallel to the plane of the sample, and fre- 
quency of AC field oj — 2irf (fixed at / = 20kHz). 
Unshunted 2D-JJAs are formed by loops of niobium is- 
lands linked through Nb—AlO x —Nb Josephson junctions 
while shunted 2D-JJAs have a molybdenum shunt resis- 
tor (with R s h ~ 2.20) short-circuiting each junction (see 
Fig.l). Both shunted and unshunted samples have rect- 
angular geometry and consist of 100 x 150 tunnel junc- 
tions. The unit cell for both types of arrays has square ge- 
ometry with lattice spacing a ~ 46/im and a single junc- 
tion area of 5 x 5/xm 2 . The critical current density for the 
junctions forming the arrays is about 600^4/cm 2 at A.2K. 
Besides, for the unshunted samples /3c(4.2if) ~ 30 and 
(3 L (A.2K) ~ 30, while for shunted samples (3 C {^-2K) ~ 1 
and (3 L (4.2K) ~ 30 where [19] (3 L {T) = 2w %° (T) and 

^(y) = 2*CjR^I c (T) ^ Here ^ Cj _ Q 58pF ig the ca _ 

pacitance, Rj ~ 10. 4£! the quasi-particle resistance (of 
unshunted array), and Ic{^-2K) ~ 150/1^4 critical cur- 
rent of the Josephson junction. $q is the quantum of 
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magnetic flux. The parameter (3l is proportional to the 
number of flux quanta that can be screened by the maxi- 
mum critical current in the junctions, while the Stewart- 
McCumber parameter (3c basically reflects the quality of 
the junctions in arrays. 

It is well established that both magnetic and trans- 
port properties of any superconducting material can be 
described via a two-component response [20], the intra- 
granular (associated with the grains exhibiting bulk su- 
perconducting properties) and intergranular (associated 
with weak-link structure) contributions [21,22]. Likewise, 
artificially prepared JJAs (consisting of superconducting 
islands, arranged in a symmetrical periodic lattice and 
coupled by Joscphson junctions) will produce a similar 
response [23]. 

Since our shunted and unshuntcd samples have the 
same value of [3l an d different values of (3c, it is pos- 
sible to verify the dependence of the reentrance effect 
on the value of the Stewart-McCumber parameter. For 
the unshunted 2D-JJA (Fig. 2a) we have found that 
for an AC field lower than 50mOe (when the array is 
in the Meissner-like state) the behavior of x'(T) is quite 
similar to homogeneous superconducting samples, while 
for h ac > 50mOe (when the array is in the mixed- 
like state with practically homogeneous flux distribution) 
these samples exhibit a clear reentrant behavior of sus- 
ceptibility [13]. At the same time, the identical exper- 
iments performed on the shunted samples produced no 
evidence of any reentrance for all values of h ac (see Fig. 
2b). It is important to point out that the analysis of the 
experimentally obtained imaginary component of suscep- 
tibility x"{T) shows that for the highest AC magnetic 
field amplitudes (of about 200mOe) dissipation remains 
small. Namely, for typical values of the AC amplitude, 
h ac = lOOmOe (which corresponds to about 10 vortices 
per unit cell) the imaginary component is about 15 times 
smaller than its real counterpart. Hence contribution 
from the dissipation of vortices to the observed phenom- 
ena can be safely neglected. 

To further study this unexpected behavior we have also 
performed experiments where we measure x'(T) for dif- 
ferent values of Hdc keeping the value of h ac constant. 
The influence of DC fields on reentrance in unshunted 
samples is shown in Fig. 3. On the other hand, the 
shunted samples still show no signs of reentrance, follow- 
ing a familiar pattern of field-induced gradual diminish- 
ing of superconducting phase (very similar to a zero DC 
field flat-like behavior seen in Fig. 2b). 

To understand the influence of DC field on reentrance 
observed in unshunted arrays, it is important to empha- 
size that for our sample geometry this parallel field sup- 
presses the critical current Ic of each junction without 
introducing any detectable flux into the plaquettes of 
the array. Thus, a parallel DC magnetic field allows us 
to vary Ic independently from temperature and/or ap- 
plied perpendicular AC field. The measurements show 



(see Fig. 3) that the position of the reentrance is tuned 
by Hdc- We also observe that the value of tempera- 
ture T min (at which x'(T) has a minimum) first shifts 
towards lower temperatures as we raise Hd c (for small 
DC fields) and then bounces back (for higher values of 
Hdc)- This non-monotonic behavior is consistent with 
the weakening of Ic(T) and corresponds to Fraunhofcr- 
like dependence of the Josephson junction critical current 
on DC magnetic field applied in the plane of the junction. 
We measured Ic from transport current-voltage charac- 
teristics, at different values of Hd c at T = 4.2K and 
found that ^{T — 4.2K), obtained from the isotherm 
T = 4.2K (similar to that given in Fig. 3), shows the 
same Fraunhofcr-likc dependence on Hd c as the critical 
current Ic{Hdc) of the junctions forming the array (see 
Fig. 4). This gives further proof that only the junction 
critical current is varied in this experiment. This also 
indicates that the screening currents at low temperature 
(i.e., in the reentrant region) are proportional to the crit- 
ical currents of the junctions. In addition, this shows 
an alternative way to obtain Ic(Hd c ) dependence in big 
arrays. And finally, a sharp Fraunhofer-like pattern ob- 
served in both arrays clearly reflects a rather strong co- 
herence (with negligible distribution of critical currents 
and sizes of the individual junctions) which is based on 
highly correlated response of all single junctions forming 
the arrays, thus proving their high quality. Such a unique 
behavior of Josephson junctions in our samples provides 
a necessary justification for suggested theoretical inter- 
pretation of the obtained experimental results. Namely, 
based on the above-mentioned properties of our arrays, 
we have found that practically all the experimental re- 
sults can be explained by analyzing the dynamics of just 
a single unit cell in the array. 

III. THEORETICAL INTERPRETATION AND 
NUMERICAL SIMULATIONS 

To understand the different behavior of the AC sus- 
ceptibility observed in shunted and unshunted 2D- J J As, 
in principle one would need to analyze in detail the flux 
dynamics in these arrays. However, as we have previ- 
ously reported [13], because of the well-defined periodic 
structure of our arrays (with no visible distribution of 
junction sizes and critical currents), it is reasonable to 
expect that the experimental results obtained from the 
magnetic properties of our 2D- JJAs can be quite satis- 
factory explained by analyzing the dynamics of a single 
unit cell (plaquettc) of the array. An excellent agreement 
between a single-loop approximation and the observed 
behavior (seen through the data fits) justifies a posteri- 
ori our assumption. It is important to mention that the 
idea to use a single unit cell to qualitatively understand 
PME was first suggested by Auletta et al. [24]. They 
simulated the field-cooled DC magnetic susceptibility of 
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a single-junction loop and found a paramagnetic signal 
at low values of external magnetic field. 

In our calculations and numerical simulations, the unit 
cell is a loop containing four identical Josephson junc- 
tions and the measurements correspond to the zero-field 
cooling (ZFC) AC magnetic susceptibility. We consider 
the junctions of the single unit cell as having capacitance 
Cj, quasi-particle resistance Rj and critical current Iq- 
We have used this simple four-j unctions model to study 
the magnetic behavior of our 2D-JJA by calculating the 
AC complex magnetic susceptibility \ — x' + *x" as a 
function of T, (3c and (3l- Specifically, shunted samples 
are identified through low values of the McCumber pa- 
rameter ((3c ~ 1) while high values ((3c S> 1) indicate an 
unshunted 2D-JJA. 

If we apply an AC external field B ac (t) = [ioh ac cos wt 
normally to the 2D-JJA and a DC field B dc = ^oH dc 
parallel to the array, then the total magnetic flux <!>(£) 
threading the four-junction superconducting loop is given 
by <I>(t) = <& ext (t)+LI(t) where L is the loop inductance, 
^ext(t) = SB ac (t) + IdBdc is the flux related to the ap- 
plied magnetic field (with I x d being the size of the single 
junction area, and S ~ a 2 being the projected area of the 
loop), and the circulating current in the loop reads 



i(t) = i c (T)^ut) + + 



Here 4>i(t) is the gauge- invariant superconducting phase 
difference across the ith junction, and $o is the magnetic 
flux quantum. 

Since the inductance of each loop is L = [i a ~ 64pH 
and the critical current of each junction is Ic — 150/jA, 
for the mixed-state region (above 50mOe) we can safely 
neglect the self-field effects because in this region LI(t) 
is always smaller than & ext (t). Besides, since the length 
I and the width w of each junction in our array is 
smaller than the Josephson penetration depth Aj = 
\J ( E > o/27r/io<i/cO (where j c o is the critical current density 
of the junction, and d = 2Xl + £ is the size of the con- 
tact area with A l (T) being the London penetration depth 
of the junction and £ an insulator thickness), namely 
I ~ w ~ 5/ito and Aj ~ 20/im (using j c o ~ 600 Aj cm 2 
and Al — 39nm for Nb at T = 4.2if), we can adopt the 
small-junction approximation [19] for the gauge- invariant 
superconducting phase difference across the ith junction 
(for simplicity we assume as usual [13] that <px — <p 2 — 
4>3 =4>4 = 4>i) 



4>i(t) = (po(Hdc) 



2irB ac (t)S 



(2) 



where 4> (H dc ) = 4> (0) + 2TTii H dc dl/$ with O (O) being 
the initial phase difference. 

To properly treat the magnetic properties of the sys- 
tem, let us introduce the following Hamiltonian 



4 i 
H(t) = J^[l - cos&(*)] + ^LI 2 (t) 



(3) 



which describes the tunneling (first term) and inductive 
(second term) contributions to the total energy of a single 
plaquette. Here, J(T) = ($o/2ir)I c (T) is the Josephson 
coupling energy. 

The real part of the complex AC susceptibility is de- 
fined as 



X '(T,h ac ,H dc ) = 



dM 
dh ac 



where 



M(T,h ac ,H dc ) = -i 1 <)H 



V \dh„ 



(4) 



(5) 



is the net magnetization of the plaquette. Here V is the 
sample's volume, and < ... > denotes the time averaging 
over the period 2-k/lo, namely 



<A> 



1 f 2 



d(ujt)A(t) 



(6) 



Taking into account the well-known [25] analytical ap- 
proximation of the BCS gap parameter (valid for all tem- 



(f) pcraturcs), A(T) = A(0)tanh(2.2 



T c -T 



for the ex- 



plicit temperature dependence of the Josephson critical 
current 



I C (T)=I C (0) 



A(T) 



A(0) 



tanh 



A(T) 



2k R T 



(7) 



we successfully fitted all our data using the following set 
of parameters: </>o(0) = -| (which corresponds to the max- 
imum Josephson current within a plaquette), /?l(0) = 32, 
(3 C (0) = 32 (for unshunted array) and (3 C (0) = 1.2 (for 
shunted array) . The corresponding fits are shown by solid 
lines in Figs. 2 and 3 for the experimental values of AC 
and DC field amplitudes. 

In the mixed-state region and for low enough fre- 
quencies (this assumption is well-satisfied because in our 
case u) ^ uolr and u <§; ujlc where ujlr = R/L and 
ujlc = 1/V LC are the two characteristic frequencies of 
the problem) from Eqs.(3)-(6) we obtain the following 
approximate analytical expression for the susceptibility 
of the plaquette 



X '(T,h ac ,H dc ) ~ ~ X o(T)[(3L(T)h(b) cos 



2H d 



Ho 



(8) 



(^) - <; 



where X o(T) = irS 2 I c (T)/V4> , H = $ /(27r/x d/) ~ 
lOOe, h(b) = Jo (26) - J 2 (2b), and f 2 (b) = J (b) - 
bJi(b)~3J 2 (b) + b.J 3 (b) with b = 2i:S^ h ac /^ a and J n (x) 
being the Bessel function of the nth order. 
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Notice also that the analysis of Eq.(8) reproduces the 
observed Fraunhofer-likc behavior of the susceptibility in 
applied DC field (see Fig. 4) and the above-mentioned fine 
tuning of the reentrance effect (see also Ref.13). Indeed, 
according to Eq.(8) (and in agreement with the obser- 
vations), for small DC fields the minimum temperature 
T m i n (indicating the beginning of the reentrant transi- 
tion) varies with Hd c as follows, 1 — T. m i n /Tc — Hd c /Ho. 

To further test our interpretation and verify the influ- 
ence of the parameter (3c on the reentrance, we have also 
performed extensive numerical simulations of the four- 
junction model previously described but without a sim- 
plifying assumption about the explicit form of the phase 
difference based on Eq.(2). More precisely, we obtained 
the temperature behavior of the susceptibility by solv- 
ing the set of equations responsible for the flux dynam- 
ics within a single plaquette and based on Eq.(l) for 
the total current I(t), the equation for the total flux 
<£>(£) = $ ex t(t) + LI(t) and the flux quantization condi- 
tion for four junctions, namely 4>i(t) = ^ (n + where 

n is an integer. Both Euler and fourth-order Rungc- 
Kutta integration methods provided the same numerical 
results. In Fig. 5 we show the real component of the simu- 
lated susceptibility x(T) corresponding to the fixed value 
of f3c(T = 4.2K) = 1 (shunted samples) and different 
values of (3 L (T = A.2K) = 1, 10, 15, 20, 30, 40, 50, 60, 90, 
150 and 200. As expected, for this low value of (3c reen- 
trance is not observed for any values of (3l- On the other 
hand, Fig. 6 shows the real component of the simulated 
X(T) but now using fixed value of (3 L {T = A.2K) = 30 
and different values of (3 C {T = A.2K) = 1, 2, 5, 10, 20, 
30 and 100. This figure clearly shows that reentrance 
appears for values of (3c > 20. In both cases we used 
h ac = 70mOe. We have also simulated the curve for 
shunted ((3l — 30, (3c = 1) and unshunted ((3l = 30, 
(3c = 30) samples for different values of h ac (see Fig. 7). 
In this case the values of the parameters (3l and (3c were 
chosen from our real 2D-JJA samples. Again, our simu- 
lations confirm that dynamic reentrance does not occur 
for low values of (3c, independently of the values of (3l 
and h ac . 

The following comment is in order regarding some ir- 
regularities visibly seen in Figs.(5)-(7) around the tran- 
sition regions from non-reentrant to reentrant behavior. 
It is important to emphasize that the above irregulari- 
ties are just artifacts of the numerical simulations due 
to the conventional slow-converging real-time reiteration 
procedures [13]. They neither correspond to any exper- 
imentally observed behavior (within the accuracy of the 
measurements technique and data acquisition), nor they 
reflect any irregular features of the considered here theo- 
retical model (which predicts a smooth temperature de- 
pendence seen through the data fits). As usual, to avoid 
this kind of artificial (non-physical) discontinuity, more 
powerful computers are needed. 



Based on the above extensive numerical simulations, 
a resulting phase diagram (3c — (3l (taken for T = IK, 
h ac = 70mOe, and H dc = 0) is depicted in Fig. 8 which 
clearly demarcates the border between the reentrant 
(white area) and non-reentrant (shaded area) behavior 
in the arrays for different values of (3l{T) and (3c (T) pa- 
rameters at given temperature. In other words, if (3l and 
(3c parameters of any realistic array have the values in- 
side the white area, this array will exhibit a reentrant 
behavior. In addition, this diagram shows that one can 
prepare a reentrance exhibiting array by changing one of 
the parameters (usually, it is much easier to change (3c 
by tuning the shunt resistance rather than the geometry 
related inductance parameter (3l). 

It is instructive to mention that a hyperbolic-like char- 
acter of (3l vs (3c law (seen in Fig. 8) is virtually present in 
the approximate analytical expression for the susceptibil- 
ity of the plaquette given by Eq.(8) (notice however that 
this expression can not be used to produce any quantita- 
tive prediction because the neglected in Eq. (8) frequency- 
related terms depend on (3l and (3c parameters as well). 
A qualitative behavior of the envelope of the phase dia- 
gram (depicted in Fig. 8) with DC magnetic field Hd c (for 
T = IK and h ac = 70mOe), obtained using Eq.(8), is 
shown in Fig. 9. 

And finally, to understand how small values of (3c 
parameter affect the flux dynamics in shunted arrays, 
we have analyzed the &tot(&ext) diagram. Similarly to 
those results previously obtained from unshunted sam- 
ples [13], for a shunted sample at fixed temperature this 
curve is also very hysteretic (see Fig. 10). In both cases, 
(&tot vs. <fr ex t shows multiple branches intersecting the 
line <&t t = which corresponds to diamagnetic states. 
For all the other branches, the intersection with the line 
&tot — &ext corresponds to the boundary between dia- 
magnetic states (negative values of \') and paramagnetic 
states (positive values of x')- As we nave reported be- 
fore [13], for unshunted 2D-JJA at temperatures below 
7.QK the appearance of the first and third branches adds 
a paramagnetic contribution to the average value of 
When (3c is small (shunted arrays), the analysis of these 
curves shows that there is no reentrance at low temper- 
atures because in this case the second branch appears 
to be energetically stable, giving an extra diamagnetic 
contribution which overwhelms the paramagnetic contri- 
bution from subsequent branches. In other words, for low 
enough values of (3c (when the samples are ZFC and then 
measured at small values of the magnetic field), most of 
the loops will be in the diamagnetic states, and no para- 
magnetic response is registered. As a result, the flux 
quanta cannot get trapped into the loops even by the 
following field-cooling process in small values of the mag- 
netic field. In this case the superconducting phases and 
the junctions will have the same diamagnetic response 
and the resulting measured value of the magnetic sus- 
ceptibility will be negative (i.e., diamagnetic) as well. On 
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the other hand, when [3c is large enough (unshunted ar- 
rays), the second branch becomes energetically unstable, 
and the average response of the sample at low tempera- 
tures is paramagnetic (Cf. Fig. 7 from Ref. [13]). 

In conclusion, our experimental and theoretical re- 
sults have demonstrated that the reentrance phenomenon 
(and concomitant PME) in artificially prepared Joseph- 
son Junction Arrays is related to the damping effects 
associated with the Stewart-McCumber parameter (3c- 
Namely, reentrant behavior of AC susceptibility takes 
place in the underdamped (unshunted) array (with large 
enough value of (3c) and totally disappears in over- 
damped (shunted) arrays. 
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FIG. 1. Left: photograph of the unshunted array; right: 
photograph of the shunted array. 
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C FIG. 4. The critical current Ic (open squares) and the real 

part of AC susceptibility x' (solid triangles) as a function of 
FIG. 2. Experimental results for X '(T,h ac ,H dc ): (a) un- DC field H dc for T = 4.2K (from Ref.13). 
shunted 2D-JJA for h ac — 10 and lOOmOe; (b) shunted 
2D-JJA for h ac . = 10, 25, and 200mOe. In all these experi- 
ments Hdc = 0. Solid lines are the best fits (see text). 
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FIG. 5. Numerical simulation results for h ac = 70mOe, 
H dc = 0, /3 C {T = 4.2*0 = 1 and for different values of 
L {T = 4.2K) based on Eqs.(4)-(7). 




^ -0.05H 



4 5 6 7 8 
Temperature (K) 



FIG. 6. Numerical simulation results for hac = 70mOe, 
H dc = 0, /3 L (T = 4.2*0 = 30 and for different values of 
f3 c (T = 4.2*0 based on Eqs.(4)-(7). 
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FIG. 7. Curves of the simulated susceptibility (for Hdc = 
and for different values of h a c) corresponding to (a) unshunted 
2D-JJA with (3 L (T = 4.2*0 = 30 and Pc(T = 4.2*0 = 30 ; 
(b) shunted 2D-JJA with (3 L (T = 4.2K) = 30 and 
(3c{T = 4.2*T) = 1. 
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FIG. 8. Numerically obtained phase diagram (taken for 
T = IK, h ac = 70raOe, and Hdc = 0) which shows the 
border between the reentrant (white area) and non-reentrant 
(shaded area) behavior in the arrays for different values of (5l 
and fic parameters. 
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FIG. 10. Numerical simulation results showing <I>t ot vs. 
$> ext for shunted 2D-JJA with f3 L (T = 4.2K) = 30 and 
Pc(T = 4.2K) = 1. 
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FIG. 9. A qualitative behavior of the envelope of the phase 
diagram (shown in previous figure) with DC magnetic field 
H dc (for T=1K and h ac = 70raOe) obtained from Eq.(8). 
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